Two-body anticorrelation in a harmonically trapped ideal Bose gas 
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We predict the existence of a dip below unity in the second-order coherence function of a partially 
condensed ideal Bose gas in harmonic confinement, signaling the anticorrelation of density fluctu- 
ations in the sample. The dip in the second-order coherence function is revealed in a canonical- 
ensemble calculation, corresponding to a system with fixed total number of particles. In a grand- 
canonical ensemble description, this dip is obscured by the occupation-number fluctuation catas- 
trophe of the ideal Bose gas. The anticorrelation is most pronounced in highly anisotropic trap 
geometries containing small particle numbers. We explain the fundamental physical mechanism 
which underlies this phenomenon, and its relevance to experiments on interacting Bose gases. 

PACS numbers: 05.30.Jp, 03.75.Hh 



I. INTRODUCTION 

The "grand-canonical fluctuation catastrophe" [1 of 
the ideal Bose gas has been the subject of longstand- 
ing attention since the early days of quantum statistical 
mechanics PH7] (for a comprehensive review see Ref . |i8j ) . 
The epithet catastrophe refers to the pathologically large 
occupation-number fluctuations predicted for the con- 
densate mode below the condensation temperature Tc 
in the grand canonical ensemble. Such a prediction is 
pathological because the fluctuations must in fact be- 
come small and eventually vanish as atoms condense into 
the ground state as the temperature T approaches zero. 
This anomaly is resolved by the introduction of weak in- 
terparticle interactions, but for a strictly ideal Bose gas 
the problem must be considered within the canonical or 
microcanonical formalisms [TJ [SHU] , corresponding to a 
diffusively isolated system with a fixed total number of 
particles, in order to obtain physically sensible results. 

The advent of Bose-Einstein condensates (BECs) in 
trapped atomic gas experiments |12f[Ti] has lead to a 
revival of interest in the condensate number fluctuation 
problem [TJ [51-fTTl [T51 [T5] , as the isolation of these sys- 
tems by the trapping potential makes them practically 
ideal realizations of the restricted statistical ensembles 
with flnite and well-defined total particle numbers. Par- 
allel developments of experimental techniques to probe 
these systems via the measurement of atom-field correla- 
tions |17ff32] such as the two-point, two-body coherence 
function |33| 

,(2) . = {^\v)i'^{v'mv'mr)) 

(*t(r)5'(r))(*t(r')^'(r'))' 

where ^^(r) and *I'(r) are the bosonic creation and an- 
nihilation field operators, respectively, allow us to refor- 
mulate the problem in terms of the spatial coherence of 
the Bose field. 



Within the grand-canonical description of the ideal 
Bose gas one finds that the second-order coherence at 
zero spatial separation is given by ^'■^■'(r, r) = 2 [33, ,34j , 
even in the limit T — > 0. This is simply the matter- 
wave analog of the celebrated Hanbury Brown and 
Twiss (HBT) effect of photon "bunching" [35, 3^ — 
a two-particle interference effect due to the fundamen- 
tal symmetry of many-body Bose wavefunctions j37j. 
This prediction, while being correct for high-temperature 
gases |19l 12^ E5] , is at variance with the expectation 
that a Bose system at temperatures well below Tc should 
exhibit higher-order coherence and thus g^'^\Y,Y) « 1 
(see, e.g., Refs. [38J |39]). This discrepancy arises di- 
rectly from the spuriously large fluctuations of the con- 
densate occupation in the grand-canonical ensemble, and 
is therefore resolved in a noninteracting gas by adopting 
a canonical-ensemble treatment of the atom- number fluc- 
tuations in the ground state [33] . or by the introduction 
of weak interparticle interactions. 

In this article we show that the canonical ensemble 
treatment of the ground state (condensate mode) occu- 
pation fluctuations also reveals an additional and coun- 
terintuitive effect: an anticorrelation dip g^'^^{v,v') < 1 
in the two-body coherence function evaluated at partic- 
ular flnitely separated (|r — r'| > 0) points r and r'. 
This is a somewhat surprising result for noninteracting, 
bosonic atoms, as density anticorrelation effects are usu- 
ally associated with either strong repulsive interactions 
in bosonic systems 001 or Pauli blocking in fermionic 
systems [22_. We show here that the anticorrelation of 
density fluctuations in the trapped ideal Bose gas arises 
from the intrinsic amplitude anticorrelation of excited 
single-particle trap eigenstates through their interference 
with the condensate mode. This effect becomes pro- 
nounced when the lowest lying excited trap eigenstates 
are highly populated, but requires also that the conden- 
sate mode itself exhibits partial second-order coherence, 
and is thus obscured by the fluctuation catastrophe in a 
grand-canonical calculation of g'^^'(r,r'). 
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The two-body coherence function of a harmonically 
trapped ideal Bose gas has previously been analyzed in 
Refs. |33l [M] , In these works, the canonical ensemble 
treatment of the condensate-mode population fluctua- 
tions was implemented only to the leading order — an ap- 
proximation that is well justified for systems containing a 
very large number {N ~ 10^) of atoms. While this is suf- 
ficient for describing the primary effect of condensation 
on g^'^\r,r') — the establishment of local second-order 
coherence ^^^•'(r, r) w 1 — the restriction to very large N 
also implies a vanishingly small anticorrelation dip. As 
we show here, the anticorrelation dip reaches appreciable 
sizes only in systems containing ~ 100-1000 atoms. In 
this regime, however, the leading-order approximation of 
the condensate population fluctuations becomes inade- 
quate for quantitatively predicting the magnitude of the 
dip. Here we treat the condensate-mode fluctuations ex- 
actly, using a numerical implementation of a recurrence 
relation for the canonical partition function [11, 41-43 
and its relation to the full counting statistics of the con- 
densate mode population [SI dni E]. Our treatment 
therefore provides a more complete and quantitatively 
accurate description of the two-body coherence function 
of a harmonically trapped ideal Bose gas. As we dis- 
cuss in Sec. |III[ the ideal-gas regime in which a canonical 
treatment of fluctuations is necessitated, and in which the 
anticorrelation dip is expected to occur, is within reach of 
modern experimental atomic-physics techniques. More- 
over, the ideal-gas analysis we present here serves as an 
instructive prototype for closely related anticorrelation 
effects in interacting-gas experiments |26| . which will be 
discussed elsewhere 1441. 



II. TWO-BODY CORRELATION OF AN IDEAL 
BOSE GAS IN A HARMONIC TRAP 

We focus here on a system of N noninteracting Bose 
atoms confined to a three-dimensional (3D) harmonic 
trapping potential with frequencies uii {i = x,y,z). We 
denote the eigenstates of the trap by Cn(r) {n > 0), and 
their corresponding energies by En- The annihilation op- 
erator Qn = J dr (r)vE'(r) destroys an atom in the state 
C„(r). In Appendix [A| we derive an approximate form 



.(2) 



(r,r') = l 



\G'^'Hr,r')\ 
p(r)p(r') 
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p(r)p(r') 



(2) 

for the second-order coherence function of the gas, ne- 
glecting higher-order terms that have little effect on 
the anticorrelation dip we study in this article. Here, 
G'^-'(r, r') = (\I'^(r)4'(r')) is the first-order coherence 
function of the gas, p{r) — G'^^^(r, r) is the atomic den- 
Po{r) — (aQao)|Co(r)P is the density of atoms in 



sity, 

the ground state, and g. 



a^alanan) / (KO'n)'^ is the 
second-order coherence of the mode Cn(r). 

All expectation values in Eq. ([2| refer, in principle, to 
traces with respect to the density matrix appropriate to 



the canonical ensemble. Our analysis is, however, sub- 
stantially simplified by calculating all occupation num- 
bers in the grand-canonical ensemble with mean total 
atom number N. Although the grand-canonical ensem- 
ble greatly overestimates the fluctuations of the conden- 
sate population, it predicts the mean occupations of the 
trap modes with very little error, compared to the exact 
canonical-ensemble values (see, e.g., Refs. [5J [10]). Fol- 
lowing Ref . |32| ) for a given atom number A'^ and temper- 
ature T, we calculate occupation numbers in the grand- 
canonical ensemble with chemical potential /x determined 
by the implicit equation 



N 



■' 1 = 1 i ^ " 



(3) 



where we have introduced the notation Ci = hiOi/kBT 
for compactness. We can therefore use the known grand- 
canonical expressions for the first-order coherence func- 
tion IMIISI 



GW(r,r') = E' 



n 



exp 



tanh(^)(^ 
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and the atomic density 
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tanh(i|^)^ 
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, (5) 



where ai = \J h/mLUi is the characteristic length of the 
harmonic-oscillator potential in the ith direction. More- 
over, the ground-state (condensate) occupation is given 
in this approach by No = {a\ao) = (e"'^/'^^^ — 1)^"^, and 
thus the ground state density 



2 / 2 



(6) 



Implicit in Eq. Q is the assumption that the number 
fiuctuations of excited modes are well-approximated by 
the usual grand-canonical predictions, and in particular 
that no correlation exists between the fiuctuations of any 
two distinct modes. The validity of this assumption for 
the purposes of exhibiting the anticorrelation dip is dis- 
cussed in Appendix [A] Equation ^ does, however, ex- 
plicitly account for the canonical-ensemble number fluc- 
tuations of the ground state, which we quantify by |47| 



(2) _ (qoQoQoQo) _ -, J_ 



(AA^o) 



(7) 



where ANq = [{Nq) — {Nq}^]^^^ is the root-mean-square 
(rms) dispersion of Nq. 
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FIG. 1: Second-order coherence Qq"' [Eq. ([t])] of the conden- 
sate mode in an isotropic 3D harmonic trapping potential 
{uJi = ui, i = x,y,z), as a function of the reduced tempera- 
ture T/T^, where Tg = T^{N) is the thermodynamic-hmit 
value for the ideal-gas condensation temperature [Eq. ||8|]. 
The three curves correspond to three distinct values of the 
total atom number A''. 



Neglecting the last term on the right-hand side of 
Eq. Q we regain the usual grand-canonical form for 
g^'^^{r,r'), from which we find, for any T > 0, a mono- 



tonically decreasing profile of the two-body coherence 
g*-^-* (r, r -I- Ar) as a function of the relative separation 
Ar — r' — r, with the coherence function decreasing 
from its peak value of g^^\r,r) = 2 at zero separation 
to g^^\r,r + Ar) = 1 at infinitely large separations |Ar|. 
While this behavior is valid above the transition temper- 
ature Tc, the prediction that (7'^^(r,r) = 2 even at tem- 
peratures T < Tc is in contrast to the partially coherent 
behavior g^^'' (r, r) < 2 expected in the presence of a Bose 
condensate [37j — a manifestation of the grand-canonical 
fiuctuation catastrophe in terms of the two-body coher- 
ence function. 

The approximate form ^ of the second-order coher- 
ence function generalizes the one employed in Refs. |33l 
that itself includes a leading-order description of the 
suppression of condensate number fiuctuations below the 
critical temperature. Specifically, the form of g(^)(r,r') 
used in those references is obtained from Eq. ^ in the 
limit that the condensate mode exhibits Poissonian num- 

(2) 

ber fiuctuations (i.e., that g^ = 1). Such an approach 
neglects the fact that the fiuctuations of the condensate 
occupation vary smoothly (in a finite system) from those 

(2) 

of a thermal mode (iJq — 2) to Poissonian as the system 
temperature is lowered through the Bose-condensation 
transition. The regime of intermediate fiuctuation statis- 

(2) 

tics 1 < 5'q < 2 becomes increasingly narrow as the to- 
tal atom number TV increases, and approaches a sudden 
step-change at the critical temperature in the thermody- 
namic limit. The approximation g^ = 1 is thus well jus- 
tified for large atom numbers, as the second-order coher- 
ence of the condensate mode deviates significantly from 
this value only over a very narrow temperature range on 
the condensed side of the transition, whereas above Tc 
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FIG. 2: (Color online) (a) Symmetrically evaluated nonlocal 
coherence function p'^' (a;, — a:) = gr'^'(a;x, — a;x) for an ideal 
Bose gas of N = 100 atoms in an isotropic 3D harmonic trap- 
ping potential {uji = uj, i = x,y,z). The three curves cor- 
respond to three distinct reduced temperatures T/Tq, where 
Tc = Tc{N) is the thermodynamic limit value for the ideal- 
gas condensation temperature [Eq. Ql. (b) The minimal 
value of g''^\x, -x) as a function of T/T^, for TV = 100, 200, 
and 400 atoms, (c) The value of the local two-body coherence 
function at the trap center, (;''^'(0,0). The upper horizontal 
axis in (b) shows T as a proportion of the finite-size corrected 
critical temperature Tc [Eq. ||9|| for = 100 atoms. The in- 
set shows the behavior of the second-order coherence g^-^^ (0, 0) 
at the trap center as T — > 0, where we find 3^2' (0,0) ^ 1-1/TV 
(see text). 



(where w 2) the spurious contribution of the last 
term in Eq. pi is neghgible anyway, as po{r) <^ pir) in 
this regime |48] . 

For smaller total atom numbers, however, the crossover 

(2) 

region becomes broad and precise knowledge of g^ is 
required to accurately calculate the two-body correlation 
function g(^^(r,r'). Accordingly, it is crucial to evalu- 
ate the variance (AA^'o)^ and hence ^q^-* with as few ap- 
proximations as possible. Here we employ a recurrence 
relation for the canonical partition function [H [101 [H] 
to numerically calculate the canonical-ensemble values of 
(AiVo)2 (and thus g[f'') exactly |49j- 
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A. Isotropic trap 

We consider first the case of an isotropic 3D trapping 
potential {ui = oj, i — x,y,z), and show in Fig. n] the 

(2) . — 

calculated values of as functions of the temperature 
T, for N = 100, 200, and 400 atoms. For these small 

(2) (2) 

atom numbers the crossover from gg « 2 to w 1 
occurs over a broad range of temperatures. We quote the 
temperature in Fig.[T]as a fraction of the thermodynamic- 
limit result for the transition temperature of the ideal 
Bose gas with N atoms 

where the Riemann zeta function C(3) ~ 1.202. The 
temperature corresponds to the onset of saturation of 
the excited-state population and is easily derived in the 
semiclassical limit fc^T ^ hui (see, e.g., Ref. [50J. In 
systems with a relatively small N, finite-size effects result 
in a lowering of the effective condensation temperature 

to 

^ _^oA 0.7275 E.^. A 

and we find Tc/T° w 0.840, 0.877, and 0.901, for N = 
100, 200, and 400 atoms, respectively. 

We note that the forms of ^q^"* presented in Fig. [l]are 
reminiscent of the behavior of condensate fiuctuations 
near the transition to condensation in the interacting 
Bose gas (see Ref. and references therein). The fact 
that such behavior is also obtained for the canonical ideal 
Bose gas is consistent with the arguments of Wilkens and 
Weiss [10 that the constraint of fixed particle number N 
in the canonical ensemble acts as an effective interparticle 
interaction that bestows a second-order character upon 
the ideal-gas transition [52'. 

In Fig.|2ja) we plot the symmetrically evaluated nonlo- 
cal coherence function g^^^ {x, —x) = g^^-' (xx, — a;x), cal- 
culated using Ens. (|2|-0 for the same parameters as the 
results of Fig. Ill An anticorrelation dip ^^^•'(x, — a;) < 1 
at finite x is clearly visible for all three temperatures 
considered. This dip can easily be understood in terms 
of two-particle interference between excited modes of the 
trap and the condensate mode. Any given excited mode 
of the trap Cn(r), n > has a spatial structure, varying 
from positive to negative with r (for simplicity we con- 
sider the real Hermite modes of the trap). Thus when 
this mode interferes constructively with the condensate 
(which has a uniform phase across its extent) at, e.g., 
a point where Cn(r) > 0, it must correspondingly inter- 
fere destructively at a point where Cn(r) < 0. This in- 
terference produces a positive density fiuctuation at the 
former point, and a negative density fiuctuation at the 
latter point, and provides the basic mechanism for the 
anticorrelation of density fiuctuations between the two 
points. However, the condensate mode must be at least 



partially coherent (^g < 2) in order for this tendency 
towards anticorrelation to overcome the (positive) corre- 
lation between fiuctuations at the two points due to the 
HBT effect. Thus the anticorrelation dip is not observed 
in a grand-canonical treatment, due to the fluctuation 
catastrophe. A more detailed discussion of this funda- 
mental mechanism is given in Appendix [B] 

In Fig. |2|b) we plot the minimum value of g^'^\x, —x) 
as a function of the temperature of the gas. For com- 
parison, the second-order coherence g'^^^ (0, 0) at the trap 
center is shown in Fig. |2]^c). By the arguments in Ap- 
pendix |B] we expect that the magnitude of the dip in 
5*^^^ (x, —x) below the uncorrelated level of unity is, at 
leading order, approximately proportional to the peak 
fractional occupation Ni/N of the first excited state 
(along x), which can be estimated to scale roughly as 
0{1/Vn) [45., and the results of Fig. ^h) are consis- 
tent with this scaling. The anticorrelation is thus more 
pronounced for systems containing smaller numbers of 
atoms, and becomes essentially negligible for 10^ 
atoms. We also note that the second-order coherence 
g'-^' (0, 0) at the trap center does not precisely approach 
unity as T — ^ [inset to Fig. |2jc)|, but in fact tends 
towards the g^'^'>{r,r) = 1 — 1/iV result appropriate to a 
Fock state of N atoms |37] , illustrating the precision to 
which we obtain the second-order coherence in this limit. 

A second dependence of the anticorrelation dip on the 
total atom number N is that the peak anticorrelation ef- 
fect occurs at progressively higher temperatures as the 
atom number is increased. This can also be understood 
on the basis that the magnitude of the anticorrelation 
dip should be approximately proportional to Ni: The 
peak value of A^i in the ideal gas occurs at increasingly 
higher temperatures approaching Tq as the atom number 
N —i' oo |45], and indeed the temperature at which the 
maximal anticorrelation appears is reasonably well pre- 
dicted by the temperature at which the product N^Ni 
(not shown) attains its peak value, though is somewhat 
skewed to lower temperatures due to the dependence of 
gg^' on T. We note also that the dip of g^^^ (x, — x) below 
unity broadens and becomes centered around increasing 
X as the temperature of the system is lowered. This is 
easily understood in terms of decreasing (positive) contri- 
butions to g(^)(x, — x) due to HBT bunching of atoms in 
excited modes, which obscure the anticorrelation effect, 
as the temperature is reduced. 

B. Highly anisotropic trap and ID systems 

We now analyze the two-body correlation function of 
an ideal Bose gas in a highly anisotropic (cigar-shaped) 
trapping potential. This is a particularly important 
case as the anisotropic confinement, for temperatures 
well below the excitation energy in the tightly con- 
fined transverse dimensions {kgT <^ tiio±, where we 
have assumed LUy — uj^ = 3> uj^), gives access to 
one-dimensional (ID) physics which can be quite dis- 
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FIG. 3: (Color online) (a) Symmetrically evaluated nonlo- 
cal coherence function g^'^\x,—x) = g^'^\xx,—xx) for an 
ideal Bose gas of A'' = 100 atoms in a highly anisotropic 
trapping potential = 127). The three curves cor- 

respond to three distinct reduced temperatures T/Tq. (b) 
The minimal value of g^^' [x, —x) as a function of T/Tg, for 
= 100, 200, and 400 atoms, in traps with aspect ratios 
u)±/ujx = 127, 223, and 395, respectively, (c) The value of 
the two-body correlation function at the trap center 5r''^'(0, 0) 
as a function of T/T^, for N = 100, 200, and 400 atoms. The 
upper horizontal axis in (b) shows T in units of the ID critical 
temperature T^^'^ for the case of A'^ = 100 atoms. 



tinct from 3D physics |271 [551 155H55] . For example, in 
such a system in the ideal gas regime, provided that 
7V[C(3)]1/2[1j^(-2^)]-3/2 ^ uj^/w^, the role of the 3D 

critical transition temperature Tc (which can be ap- 
proximated to leading order by T^) is reduced to sig- 
naling the onset of the so-called transverse condensa- 
tion |57J. In this regime, atoms accumulate in the trans- 
verse ground state below Tc due to the saturation of pop- 
ulation in the transverse excited states, yet the system 
remains noncondensed with respect to the longitudinal 
states until the effective ID condensation temperature 



(ID) 



Nhuj^/[kBH'2-N)] 051157] is reached. 



Examples of the two-body correlation function 
g^^-* (x, —x) symmetrically evaluated along the long axis 
of a highly anisotropic trap are shown in Fig. |3ja). In all 
three cases shown we choose uj±_ to maintain the con- 
dition ksT^^'' /hLu± = 0.2, so that the system is in 
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FIG. 4: Second-order coherence Qq of the condensate mode 
as a function of T/Tq, for the same parameters as 



the quasi-lD regime at temperatures near r^^°^ and be- 
low (we find T^^^V^c « 0.230, 0.220, and 0.212 for 
N = 100, 200, and 400, respectively). 

We observe [Fig. [sj^b)] that the strongest anticorrela- 
tion occurs near T ~ T^^\ in contrast to the isotropic 
case in which it occurs near T ^ Tc- For the same 
total number of atoms N, the maximal magnitude of 
the anticorrelation dip is larger in the highly anisotropic 
trap than in a spherically symmetric trap, and as in the 
isotropic case the dip is deeper for smaller N. In Fig. |4] 
we show the values of the second-order coheience func- 
tion gl^^ used to calculate the results of Fig. bl 

J2) . ^ , . (2) 



We note 

that the crossover from ^q^^ ss 2 to g^"^ « 1 occurs over 
a very broad range of temperatures, indicating the im- 

(2) 

portance of using precise values for g^ to accurately 
evaluate (7'^^^(r,r') in this quasi-lD regime. The result- 
ing behavior of the second-order coherence i?*^^-' (0, 0) at 
the trap center is shown in Fig. |3jc). 

In Fig. ^ we plot the full matrix structure g^^'(a;,a;') = 
g'^^^ (xx, arx) of the coherence function evaluated on the 
long axis of the anisotropic trapping potential, for the 
case = 100. Along the diagonal x' = x we observe 
the suppression of g^-^^x^x) from the HBT value of 2 
(white) at large \x\ down to ^'■^■'(O, 0) « 1.5 at the origin, 
due to the presence of the condensate. The anticorrela- 
tion lobes where g^^''{x,x') < 1 (dark shades of red/gray) 
are clearly visible, and are centered on the anti-diagonal 
x' — —X at « a^, as expected (see Appendix [B| . 
We note, however, that some (comparatively small) an- 
ticorrelation can also be seen at, e.g., {x,x') = {0,ax)i 
which is evidently due to the intrinsic amplitude anti- 
correlation of even (excited) Hermite modes, illustrating 
that it is not solely the first Hermite (dipole) mode which 
contributes to the density anticorrelation feature. 

Finally in this section, we note that we have ap- 
plied our methodology to the case of iV = 100 atoms 
in a purely ID trap, and compared the results for 
min{(/(^' (x, —x)} (which agree very closely with the plot- 
ted case of u!±/ujx — 127) with the correlations of this sys- 
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FIG. 5: (Color online) Structure of the two-body coherence 
function g'^' (a;, a;') = g^^\xx., x'x.) on the long axis of a highly 



anisotropic trap. Parameters are A'^ 
127, and T/T^ = 0.2 (cf. Fig.^. 



100 atoms, U^/uJa: 



tern estimated using Monte-Carlo sampling of the canon- 
ical ensemble in the coherent-state representation [55|. 
These numerical results show qualitative agreement with 
our semi-analytic calculations, and in fact suggest that 
the maximal suppression of g'^^\x, —x) is ~ 20% larger at 
its peak (as a function of T), and that the anticorrelation 
subsides somewhat more slowly with increasing temper- 
ature, due to higher-order contributions to g^'^'>{x,x') in 
the canonical ensemble (see Appendix [A| . 



III. EXPERIMENTAL CONSIDERATIONS 

The density anticorrelation we discuss in this article 
is an intrinsically in-trap effect and could therefore be 
detected, in principle, via an in situ measurement of 
(/^^^(r,r') (as proposed, e.g., in Ref. [10] )• Moreover, as 
the momentum-space correlations of the trapped ideal 
gas can be obtained (up to a normalization factor) from 
its position-space correlations by the simple substitution 
ri — {—!^)ki |46| . where and ki are the ith compo- 
nents of the position and momentum vectors r and k, 
respectively, the effect can equally well be regarded as 
an in-trap momentum-space density anticorrelation. As 
is well known, position correlations measured in far field 
after releasing a trapped gas from its confinement reflect 
(in the limit that interactions can be neglected during the 
expansion) the corresponding in situ momentum correla- 
tions of the original trapped sample [SQ- In feet, in the 
harmonically trapped non-interacting systems we have 
considered here, the position-space correlations of the gas 
are, at all times during the expansion, given by a simple 
rescaling of the in-trap position-space correlations ^34j. 



As such, the density anticorrelations we predict for the 
ideal gas should be observable in near-field expansion as 
well. This particularly simple behavior is of course pe- 
culiar to the case of an ideal gas, so we now briefly con- 
sider the prospects for realizing the ideal-gas limit exper- 
imentally, focusing on the quasi-lD regime (Sec. II B I, in 
which the anticorrelation dip is most pronounced. 

In a weakly interacting Bose gas in a strongly 
anisotropic (elongated) trapping geometry, the finite- 
size condensation behavior |45j comes into competition 
with an interaction-driven crossover to a quasicondensate 
regime, characterized by suppressed density fiuctuations 
but lacking significant nonlocal phase coherence, with a 
characteristic crossover temperature |60j 



y(lD) 



SNhuJx 



ks \'n.{Nh?(jJx/mg'^ 



(10) 



where g sa 2hLo±a is the effective ID coupling con- 
stant [HI] with a the 3D s-wave scattering length. To 
access the effectively ideal-gas regime, the system pa- 
rameters must be such that Tco^'^ ^ T^^\ in which 
case the system will be well described by the ideal- 
gas treatment of this article in the temperature range 



(ID) 



<T <T, 



(ID) 



As shown in Ref. [60 , to access this 



regime we require a < 0.01 nm, which is far smaller than 
the natural scattering lengths of the alkali-metal atoms; 
however, this extremely weak scattering could be realized 
by exploiting Feshbach-resonance techniques |62| . 

The observation that density anticorrelations can arise 
even in an ideal Bose gas is particularly relevant in light 
of recent theoretical and experimental developments con- 
cerning the observation of density correlations in time-of- 
flight, following the release of a low-dimensional (quasi- 
condensed) gas from its conflnement |23| [63] (see also 
Ref. t64i)- The central observation of these works is 
that the (axial) coherence g^'^^x.x') of the expanding 
gas develops "ripples" at an intermediate stage of the ex- 
pansion, which correspond directly to the in situ phase 
fluctuations of the original sample. As the expansion 
continues, g^^\x,x') approaches its far-field limit, which 
reflects the in-trap momentum correlations of the field. 
In Ref. 123] , these momentum correlations are calculated 
within a local-density approximation, and the far-field 
limit of the position-space coherence function is there- 
fore found to decay monotonically with |x — a;'|. 

Although the ideal-gas results presented in this ar- 
ticle are not directly applicable to interacting systems, 
they demonstrate the potentially significant infiuence of 
finite-size and inhomogeneity effects on experimentally 
measured correlations. Indeed, although density fluc- 
tuations are strongly suppressed in quasi- ID interacting 
systems such as that considered in Ref. [23 , phase fluc- 
tuations remain significant in such regimes. In a har- 
monically trapped quasicondensate, these phase fluctu- 
ations exhibit a specific spatial structure resulting from 
the inhomogeneity of the quasicondensate density pro- 
file [65i |66] . Accordingly, we might expect such systems 
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to exhibit in-trap momentum-space density anticorrela- 
tions |44| . which should manifest in position-space mea- 
surements made in far field, and may potentially yield 
corrections to the time development of correlations dur- 
ing free expansion [63 . 

More recent experiments |26j have considered the 
correlations of interacting atomic bosons in both the 
strongly anisotropic quasicondensate regime, and in 
weakly prolate traps in which the physics of the gas is 
expected to correspond more closely to the compara- 
tively transparent case of strictly three-dimensional Bose- 
Einstein condensation. In particular, anticorrelation dips 
similar to those discussed in this article were observed in 
time-of-flight expansion of a gas released from 1 : 4.75 
aspect-ratio confinement, at temperatures T O.OTc. 
Although further analysis would be required to disentan- 
gle the effects of interactions and finite expansion times 
on the measured correlations, it nevertheless seems nat- 
ural to associate these experimental results with the fun- 
damental physical mechanism discussed in this article: 
amplitude anticorrelations of excited modes yield density 
anticorrelations upon interfering with the condensate. In 
particular, at temperatures close to Tc, we expect that 
the effects of interactions are relatively subdued, serving 
at leading order to simply restructure the single-particle 
excitations of the trap into Hartree-Fock-like modes |67| . 
which should indeed exhibit the same fundamental am- 
plitude anticorrelations that underlie the density anticor- 
relations we have considered here. 



IV. SUMMARY 

In summary, we have predicted the existence of two- 
body anticorrelations in an ideal Bose gas with a fixed, 
finite total atom number in harmonic confinement. The 
corresponding dip below the uncorrelated level of unity 
in the second-order coherence function of the Bose field 
reaches experimentally resolvable magnitudes in highly 
anisotropic traps containing small atom numbers {N < 
1000). The ideal-gas regime in which the effect is 
pronounced could be reached in contemporary atomic 
physics experiments by exploiting Feshbach resonance 
techniques. Moreover, the fundamental physical mech- 
anism underlying the anticorrelation should also be rele- 
vant to experiments on interacting Bose systems. 
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Appendix A: General form of the second-order 
coherence function 



The second-order coherence function is given in general 



by 



G(2)(r,r') = (*1'(r)^'^(r')*(r')*(r)) (Al) 
= X! C(r)C(r')Cp(r')C9(r)(aLanapag>, 

mnpq 

where (• • • ) = Tr{ • • • /5} denotes a trace taken with re- 
spect to the appropriate thermodynamic density ma- 
trix. In the grand-canonical ensemble the density matrix 
is of the Gaussian form pa = Zg^e^^^"^^-'/'^"^, with 

H = Y^n^nal^CLn and N = X]„"n«ni ^ud where Zq is 
the grand-canonical partition function, whereas the den- 
sity matrix appropriate to the canonical ensemble can be 
written as 

pc^Z^'e'^/^^Pi,, (A2) 
where the projector 

Pn= I^o,A^i,---)(^o,^i,---|, (A3) 

here expressed in terms of Fock states \Nq, Ni, ■ ■ ■) over 
the basis modes Cn(r), effects the constraint to a fixed 
total number of particles N, and Zc is the canonical 
partition function. As [pc,N] = [pc,N] = 0, in both 
the canonical and grand-canonical ensembles all so-called 
anomalous averages {{am), {amCin), (am^nflp), etc.) van- 
ish, and all two-body correlations of the field are there- 
fore completely specified by the one-body correlations 
(aLa«) 



mnpq 



{al^aidpaq) 



(A4) 



- {dl^ap){aiaq) - (aj„a,)(4ap) 



= {alnoldpaq) 



which quantify the deviation of the field fiuctuation 
statistics from the Gaussian limit. In the grand-canonical 
ensemble, these cumulants vanish identically, so the 
second-order coherence function is specified completely 
by the occupation numbers Nm, and can be written |35] 



G(2)(r,r') 



(^t(i.)#(r))(,j,t(i.')^r(r')) 
+ (*^(r)*(r'))(*^(r')*(r)). (A5) 



In the canonical ensemble we have, by the form of 



the density operator pc [Eq. (A2|], Cmnpq — unless 



{■m,n) — {p,q) or {m,n) — {q,p). Thus in the case 
m = rt, a nonzero cumulant is found only when p — 
q = m. Therefore, the second-order coherence function 
in the canonical ensemble [G'^-* (r, r')]^ = [G'^''(r, r')]^-!- 
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AG'(2)(r,r'), where 
AG(^)(r,r')= Yl C(r)C:(r')Cp(r')C,(r)C„ (A6) 

mnpq 

^ ^ |Cm(^)| |Cm(^ )| ^mmmm 

+ EE [\Ur)\'\Ur')\'C„ 



-'mnnm 



C(r)C(r')C™(r')Cn(r)a 



mnmn 



To proceed, we separate off terms in Eq. ( A6 ) that involve 



the condensate mode Co(r), and use the symmetries of the 
cumulants Cmnpq Etnd the relations 

(A7V™AiV„) = {{N„, - {N^)) {N,, - (7V„))) 

^rammra -^mi-^m ^" 1) ^ ^ 



(A7) 



to obtain (after some algebra) the exact form 



G(2)(r,r') 



(*t(r)*(r))(#t(r')5,(r')) + {¥ {r)i'{r')){i'Hr')^{r)) + [((ATVo)') - A^o(A^o + 1)] |Co(r)nCo(r')l' 
+ ^'(A7VoA7V„) |Co(r)Cn(r') + Co(r')C™(r)|' 

m 

+ [{(AN^f) - N^N^ + 1)] |C™(r)nC™(r')P 

m 

+ E' E'(^^™^^") [ICn(r)nC„(r')P + C(r)C(r')C™(r')Cn(r)] , (A8) 



where Yl'm denotes a sum which excludes the condensate 
mode (m = 0) . Terms on the first line of Eq. ( A8 1 con- 



stitute the approximate form for the coherence function 
that we use in the main text of this article, from which 
the normalized form [Eq. (|2|] is obtained as f;*^^^(r,r') = 
[G(2)(r,r')]c/[G(iHr.r)G(M;r',r')], where G(i)(r,r') = 
(^■'■(r)5'(r')) [33]. This approximate form for g(^)(r,r') 
generalizes those given in Ref. [33] by allowing for fluc- 
tuation statistics of the condensate occupation between 
the limiting cases of a purely thermal or purely coherent 
state. 



The terms on the second line of Eq. ( A8 1 correspond 



to the effects of explicit correlations between the fluc- 
tuations of the condensate occupation and those of the 
excited-mode populations. Politzer [5] has shown that 
these correlations are negative, i.e., that fluctuations 
of the excited-mode populations are anticorrelated with 
those of the condensate occupation. The anticorrelations 
between mode occupations are most pronounced near the 
Bose-condensation transition, and the largest such anti- 
correlation is that between the condensate and the (de- 
generate) first excited state, for which —{ANqANi) < 
O.lA^o-^i for TV = 100 (or more) atoms in an isotropic 
trap [S". The spatial term multiplying {ANoANm) in 
Eq. (A8l is manifestly non-negative, and these anticorre- 



lations therefore yield a small additional suppression of 
the second-order coherence function relative to Eq. (|2|; 
i.e., they do not subtract from the anticorrelation dip we 
consider in this article, and can thus be safely neglected. 



The terms on the third line of Eq. ( A8 ) correspond to 



r 



excited modes from their grand-canonical levels. In fact, 
in the canonical ensemble these fluctuations are some- 
what suppressed below the grand-canonical predictions 
— an effect easily understood in terms of the effective- 
interaction model of the canonical-ensemble constraint 
discussed by Wilkens and Weiss [1^ — and these terms 
therefore do not subtract from the anticorrelation dip 



either. The final line of Eq. ( A8 1 corresponds to the cor 



relations between excited-mode population fluctuations. 
These correlations are again negative, and although the 
spatial term multiplying (AA^^AiVn) is not strictly pos- 
itive, these anticorrelations are weak (smaller in magni- 
tude than (A7VoAA''i) fS]), and yield only a small quan- 



the deviations of the number fluctuations of each of the 



titative correction to '{r,r'). 



Appendix B: Physical origin of the effect 

The physical origin of the anticorrelation dip in 
the second-order coherence function is the intrin- 
sic amplitude anticorrelation of the excited energy 
eigenstates of the trap. The underlying mecha- 
nism can be illustrated using a simple two-mode 
model, consisting of the flrst two energy eigenmodes 
of a ID harmonic trap; i.e., the ground mode 
(o(a;) = {^a^)-^^'^exp{-x^/2al) and the "dipole" 
mode Ci{x) = {2/ ^aj)^/'^{x/ar) exp{~x'^ /2al). We 
make minimal assumptions as to the correlations of the 
two modes: we assume that there are no amplitude or 
number correlations between the two modes ((ajai) ~ 
and (ajaoaldi) = (aQCio) {a\ai) , respectively), and 
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that the excited mode is incoherent, {a[a[aiai) 
2(a|ai)^ = 2Nf, while allowing for partial coherence 
of the condensate mode, (aJaQcioao) = So^^Qo^o)^ = 
q'^^Nq. With these basic assumptions, one finds that 
the second-order coherence function is given (exactly) 
by Eq. §, where G(i)(x,x') = ELo 
p{x) = G^{x,x), and po{x) — NqQ{x). We consider the 
second-order coherence between the turning points 
of the dipole mode (between which the amplitude anticor- 
relation is largest): we have Co(±crj;) = (er^ j ^ixax)^!"^ , 
and C,i{a,^) = -C,i{-a,j:) = {2e-^ / ^(jr^f/'^ , and 
thus 
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(2) 



, ^^ , , {N,-2N,f-{2-g'i^)Nl 



(A^o + 2iVi)2 



Clearly, in the case that the condensate mode is in- 
coherent {g^^ = 2), g^^-* (cr^, — (Tj;) > 1. To obtain 

g^'^\(JxT—(^x) < 1 we in fact require ^q^'' < 2 — (1 — 
2Ni/NqY . Thus partial coherence of the condensate is 
required in order for the dip to arise, which explains 
why the anticorrelation is only visible for T < Tc, and 
is obscured in a grand-canonical treatment. Moreover, 
as Ni becomes small compared to Nq, the value of 
required to observe the dip becomes smaller, consistent 
with the observation that the anticorrelation effect weak- 
ens as T — > 0. Considering the case that the condensate 
is perfectly coherent ' = 1, in the limit A^i <C A^o 



we obtain g'-'^^a^, -a^) w 1 - 4iVoiVi/(7Vo + ZA^i)^ w 
1 - 4:NoNi/N^ « 1 - ANi/N (where iV = 'Z.N,) 
— which gives a useful characterization of the depen- 
dence of the magnitude of the dip on the occupations 
of the two modes. This simple model obviously gen- 
eralizes to the case of the gas in thermodynamic equi- 
librium, in which higher trap modes also produce den- 
sity anticorrelations upon interfering with the conden- 
sate — though the dipole modes are expected to yield 
the dominant contribution to the effect, due to the rel- 
ative magnitude of both their antinodes and their occu- 
pations at equilibrium. Finally, we note that the anti- 
correlation effect can easily be understood in terms of 
the classical picture of HBT interferometry of optical 
fields [W: Consider the field a{x) radiated by an in- 
coherent source a {{\a\'^) ~ 2(|ap)2), at two points xi 
and X2 between which a{x) exhibits a tt phase difference 
[0(2^2) = — ct(a;i)]. In the presence of a second field /3{x), 
uncorrelated with a{x), that exhibits zero relative phase 
between xi and X2, the measured intensities at the two 
points will exhibit an anticorrelation (I1/2) < {Ii){l2), 
provided that the field /? exhibits sufficient coherence 

gf ^ {\m/m')' < 2- - {\mr/m'r to 

overcome the HBT bunching effect. The only differ- 
ence in the present case of a trapped Bose gas is that 
no propagation is necessary to observe the effect, as the 
phase differences between spatial points arise naturally 
from the structure of harmonic-oscillator energy eigen- 
modes. 
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